Nanometer-resolved collective micromeniscus oscillations through optical diffraction 
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We study the dynamics of periodic arrays of micrometer-sized liquid-gas menisci formed at super- 
hydrophobic surfaces immersed into water. By measuring the intensity of optical diffraction peaks 
in real time we are able to resolve nanometer scale oscillations of the menisci with sub-microsecond 
time resolution. Upon driving the system with an ultrasound field at variable frequency we observe 
a pronounced resonance at a few hundred kHz, depending on the exact geometry. Modeling the 
system using the unsteady Stokes equation, we find that this low resonance frequency is caused by 
a collective mode of the acoustically coupled oscillating menisci. 

PACS numbers: 47.55.dd,43.35.+d,42.25.Fx 



Supcrhydrophobic surfaces have attracted much the- 
oretical and experimental interest [![. The entrapment 
of gas at the textured surface reduces the actual inter- 
facial area between the solid and the liquid, which is at 
the origin of the distinctive properties of these materials, 
including the large contact angle, the low contact an- 
gle hysteresis [IJ, the self-cleaning effect [1], and the large 
hydrodynamic slip[3]. The superhydrophobic state is 
closely related to the shape of the microscopic liquid-gas 
interfaces spanning between the ridges of the texture [5( . 
In this letter we characterize the dynamics of these mi- 
cromenisci and introduce an optical diffraction measure- 
ment that reveals their nanoscale motion. We measure 
in real time the optical diffraction intensity from a peri- 
odic array of micromenisci. Driving the system with an 
ultrasound field at variable frequency, we measure its fre- 
quency response and we identify a well-defined resonance 
peak with a center frequency well below the expectations 
for a single micromeniscus. Modeling the system using 
the unsteady Stokes equation and monopole interaction, 
we show that this frequency reduction is due to acoustic 
coupling between menisci. 

FigQ] shows the geometry of our system. It con- 
sists of 1 x 1mm 2 wide hexagonal and square arrays 
of micrometer-sized cylindrical holes (radii R = 2/im 
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FIG. 1: a) Schematic figure of a single micromeniscus. b) Pat- 
tern of an array of micromenisci. Arrows denote the primitive 
translations and nearest neighbor distance a. 
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or 3/xm, depth H = 15/im, nearest neighbor distance 
a = 15,um or 25/im). The samples were fabricated from 
Si (110) using standard micro lithography. Subsequently, 
the surfaces were hydrophobized by vapor deposition of a 
monolayer of lH,lH,2H,2H-Pcrfluorodccyltrichlorosilanc, 
following [(| . The advancing and receeding contact angles 
on an unstructured surface are 7 a = 116° and 7 r = 104°. 
Upon immersing the samples into demineralized water, 
ambient air is entrapped in every hole, leaving a water- 
air meniscus behind that is pinned at the ridge of the 
hole. Owing to the hydrostatic pressure, the menisci are 
bent inwards with equilibrium curvature kq = pgh/a, 
where g is the acceleration of gravity, a is the water sur- 
face tension, and h ps 0.1m is the distance between the 
sample and the free water surface. This implies that 
the system is in diffusion equilibrium and the gas pres- 
sure in the hole is the ambient pressure. An Ar-ion laser 
(A = 488nm, s-polarized) is used to illuminate the sam- 
ple under an angle typically between 60° to 70° with 
respect to normal incidence (see Figf2J). The diffracted 
intensity is measured with a photodiode positioned at a 
selected diffraction peak, typically chosen in the vicin- 
ity of the specular reflected beam. A broadband piezo- 
electric ultrasound transducer is placed at its focal dis- 
tance from the sample. The ultrasound transducer is 
excited to emit finite wavetrains using an arbitrary func- 
tion generator. The ultrasound pressure at the sample 
is of the order 10 2 to 10 3 Pa, which is small compared 
to a critical static pressure above which filling occurrs 
P c = 2cr cos(7 a ) / R ps 2.1 • 10 4 Pa. To check the dynamic 
stability of the menisci, we increased the ultrasound pres- 
sure to much larger values and observed how the inten- 
sity oscillations disappeared at a defined threshold. The 
ultrasound pressure is kept constant during a frequency 
sweep by controlling the driving voltage according to the 
transducers frequency response. Fig|2}: shows typical raw 
data corresponding to the beginning of a wavetrain. Af- 
ter a transient lasting for a few oscillation cycles, the sig- 
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FIG. 2: a) Schematic figure of the experiment, b) Inverted 
photograph of the diffraction pattern (■& ~ 60°). Numbers 
indicate Miller indices of diffraction orders, c) AC component 
of the light intensity measured in the (1,0) diffraction order 
at the beginning of an ultrasound wavetrain. 
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FIG. 3: a) Time-resolved intensity of the (1,0) diffraction 
order, corresponding to small ultrasound pressure (AP « 
400Pa), solid line, and large ultrasound pressure (AP « 
800Pa), dashed line, b) Calculated intensity of the first 
diffraction order as a function of meniscus deflection for a rect- 
angular surface profile with groove width w = 2R (i? = 66°). 
The data is displayed simultaneously in terms of the displace- 
ment £ and the non-dimensionalized curvature kR. Both are 
related geometrically by C = 2/k(1 — y/l — R 2 k 2 /4). The 
equilibrium curvature kqR = 0.04 corresponds to fn = 30nm. 



nal becomes sinusoidal with a constant amplitude. This 
amplitude is extracted from the raw data by calculating 
the root mean square. 

The diffracted intensity depends in a highly nonlinear 
way on the deflection of the menisci. To assure a lin- 
ear relation between the actual deflection and the mea- 
sured intensity, the menisci oscillations have to be small. 



This can be seen in Fig[3^,. While the diffracted inten- 
sity follows the sinusoidal driving pressure at low driving 
amplitudes (solid line), it is distorted at larger driving 
amplitudes (dashed line). Prior to an experiment, we re- 
duce the ultrasound pressure until the undistorted sinus 
is observed. 

To find the magnitude of the corresponding meniscus 
deflections, we consider the optical diffraction of the sam- 
ple. In the Fraunhofer limit, the diffracted intensity is 
proportional to the intensity scattered by a single unit 
cell, and for each unit cell, the scattered intensity is gov- 
erned by the interference of the elementary waves emitted 
from the cell volume 0] ■ In our case the incident angle is 
large and no light reaches the bottom of the holes such 
that the interference takes place between the waves emit- 
ted from the plane silicon surface and the waves emitted 
from the meniscus. Qualitatively, as the meniscus is de- 
flected - consider the position £ of the apex of the menis- 
cus - the intensity of a diffraction order (with diffraction 
angle « incident angle) changes sinusoidally with a pe- 
riod T = A/(2ncos(i?)), where n = 1.33 is the refractive 
index of water. To analyze these simple observations in 
detail, we performed a diffractive optics calculation us- 
ing the multilayer rigorous coupled wave analysis in the 
formulation of [8( . This method allows for calculating an 
exact solution to the Maxwell equations for the optical 
response of arbitrary periodic surface profiles. In Figj3j3 
we show the resulting diffraction intensity as a function 
of the meniscus deflection. The typical distance between 
two adjacent peaks corresponds to the period evaluated 
from the simple Fraunhofer arguments above. The re- 
sult shows that the diffracted intensity is indeed linear in 
the meniscus deflection in a range [— £;.. + (i] around the 
meniscus equilibrium position Co and we find Q « 90nm. 
Note that the extent of the linear range depends on the 
incident angle It is larger for larger incident angles, 
as can be seen readily from the simple expression for T. 
Thus, the large angles that are used in experiment allow 
for large meniscus oscillations. For incident angles above 
the angle of total reflection between water and air 48.6°, 
in addition the relative contribution of the menisci to the 
scattered intensity is large. 

The theoretical result for the linear range is the key to 
convert the measured intensity variations into absolute 
meniscus deflections. Since we have to assure linearity 
between intensity and meniscus deflection at all ultra- 
sound frequencies, the peak of the resonance curve shown 
in the following has the height C,i . The slope of the linear 
range together with the relative noise of the photodiode 
determines the resolution of the deflection measurement. 
It is of the order lnm under given conditions. 

Fig[4] shows a typical measured frequency response. 
The sample displays a single resonance at f r = 153 ± 
5kHz. Similar curves were obtained for all samples. Ta- 
ble U shows that the observed resonance frequency in- 
creases both with increasing nearest neighbor distance 
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FIG. 4: Frequency response of an array of micromenisci with 
hexagonal pattern, a = 15^m and R — 3/im. Crosses show 
experimental data. The dash-dotted line shows the theory for 
a single meniscus. Theoretical data for the array is displayed 
in terms of the mean (solid line) and the root mean square 
deflection (dashed line). 
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TABLE I: Resonance frequencies of the lowest collective mode 
f c as a function of lattice constant, pattern geometry and 
menisci radius. 'H' and 'S' denote hexagonal and square pat- 
tern, respectively. f s denotes the resonance frequency of a 
corresponding single meniscus, f r the experimental result. 



and with decreasing hole size. 

To understand the observations quantitatively, we con- 
sider first the response of a single meniscus under the in- 
fluence of the pressure field P(f) = P + APe l27rft , where 
/ is the ultrasound frequency, Pq is the ambient pres- 
sure and AP is the amplitude of the ultrasound pressure. 
Since Co <C AC, R, as follows with the values evaluated 
above, we approximate the meniscus as flat in equilib- 
rium and its deflections as parabolic £(r) = ((l — r 2 /R 2 ). 
The parabolic shape implies that the curvature is uniform 
up to 0(C 3 ) giving right to the Laplace law. As described 
in the smallness of the deflections A£ <C P, and the 
high frequency / ^> v/R 2 of the oscillations (v kine- 
matic viscosity), allow us to omit the non-linear term in 
the Navier-Stokes equation and hence the dynamics of 
the system is governed by the unsteady Stokes equation. 
The system can then be described as a harmonic oscilla- 



+ K} AC* = -AP* 



(1) 



The transfer functions <j> and W account for the inertia 
and the viscous damping due to the oscillatory flow fields. 
All quantities are non-dimenionalized lo* = 2TtfR 2 /v, 
A(* = C/R and AP* = APR 2 J(pv). Since polytropic 
and thermal dissipative effects [10( can be neglected, the 
potential term reduces to a dimensionless spring constant 
K = R/{pv 2 ) (P R 2 /(2H)+4a). The first term with 
the ambient pressure is due to the isothermal compres- 
sion of the gas and is negligble throughout this work, and 
the second term with the surface tension a is due to the 
surface energy of the liquid-gas interface. The compu- 
tation of 4> and \& is performed by solving the unsteady 
Stokes equation in cylindrical coordinates with classical 
no-slip and free slip boundary conditions at the solid- 
liquid and the liquid-gas interface, respectively. A de- 
tailed account of the calculations will be given elsewhere 
While analytical expressions can be found both for 
the high and for the low frequency limit, $ and \P have to 
be computed numerically in the intermediate frequency 
range 10° < lo* < 10 2 , which is relevant for the present 
experiments. The dash-dotted line in FigH] shows the so- 
lution for a R = 3/im hole with physical parameters of 
water p = 10 3 kg/m 3 , v = l(T 6 m 2 /s, a = 73 • l(T 3 N/m 
and ultrasound pressure AP = 390Pa as used in the re- 
spective experiment. The obtained resonance frequency 
is approximately five times larger than the one observed 
experimentally. Clearly, the single meniscus theory fails 
to describe the dynamics of the system. 

To resolve the discrepancy, we consider the dynamic 
coupling between the menisci. As the menisci oscil- 
late in the external pressure field, they emit pressure 
waves that affect the other menisci. Denoting the non- 
dimcnsionalizcd distance between the i-th and the j-th 
meniscus by d*j = dij/R, the additional force acting on 
the i-th meniscus can be expressed in terms of a multi- 
pole expansion -A(*io* 2 /(4:d* 3 ) + 0{d*f 2 ), where AQ* 
is the deflection amplitude of the j-th meniscus (see e.g., 
[l^ ]. [HI). To analyze the dynamics of the entire menis- 
cus array, we extend the equation of motion of the single 
meniscus, Eq.JT]), by the additional forces generated by 
all other menisci keeping only the monopole term. We 
arrive at the coupled equations of motion 



{$(w*) + + K} AC* = -AP* 



N-M ^ 2 



(2) 

The coupling term gives rise to an additional effective 
mass, which reduces the resonance frequency, as required. 
We solve Eq[5] for the individual deflection amplitudes 
AC;* by numerical matrix inversion. From the result we 
evaluate the mean (AC*) = | Y^fJi A&*|/JVM and the 
root mean square (AC* 2 } 1/2 - I A&* \ 2 )^ 2 /NM, 
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deflection amplitude (see discussion below). The results 
are shown in Fig[4]with the solid and dashed line, respec- 
tively. For the moment note that both curves are nearly 
identical up to frequencies well above the lowest reso- 
nance. The theoretically obtained resonance frequency 
f c = 159kHz is in excellent agreement with the experi- 
mental data, showing that collective effects are crucial for 
the dynamics of the system. Moreover, the good agree- 
ment shows that monopole interaction is effective beyond 
its obvious range of applicability where o>i?. 

As can be seen in the inset of FigHJ the theory overes- 
timates the amplitude of the resonance, presumably due 
to the neglectance of the bulk dissipation that arises in 
the collective flow field. Second, and more interestingly, 
the theoretical data for the root mean square displays a 
second resonance at 290kHz which is absent in the exper- 
imental data. To understand the latter effect we plot the 
calculated amplitude and phase for the two lowest reso- 
nances as a function of the meniscus position within the 
array in FigO For the lowest resonance all menisci oscil- 
late essentially in phase, whereas at the second resonance, 
the menisci in the center of the array and the ones along 
the edge oscillate 180° out of phase and the amplitude 
displays a node at the ring shaped boundary between the 
two regions. To evaluate the diffraction intensity for such 
arrays of non-identical scatterers, we note that the varia- 
tion of the meniscus deflection gives rise to phase differ- 
ences between the waves emitted from different unit cells 
that are much smaller than 2tt - owing to the particular 
experimental condition assuring that the diffracted inten- 
sity is linear in the meniscus displacement. Extending the 
above Fraunhofer picture, one shows that the diffraction 
intensity is linear in the individual menisci deflections, 
and thus the experiment measures the mean deflection. 
Thus, we have to compare the experiment to the theoreti- 
cal mean deflection, where the second resonance is indeed 
nearly invisible. Note that the expression for the mean 
as given above accounts for the phase since the deflection 
amplitudes A£* are complex. 

The reduction of the resonance frequency predicted 
by the theory is confirmed for all surface patterns in- 
vestigated in experiment (see Table |TJ . Since the cou- 
pling between the menisci is inversely proportional to 
the distance, the frequency reduction is more pronounced 
for smaller lattice constants. Similarly, it is more pro- 
nounced for the hexagonal lattices than for the square 
ones, owing to the larger number of nearest neighbors. 
At this moment we have no clear explanation for the 
slight overestimation of the resonance frequency for the 
samples with smaller radius. 

In conclusion, highly mobile micromenisci are present 
at the textures of super hydrophobic surfaces. Their dy- 
namics are determined by collective modes with reso- 
nance frequencies that are much smaller than the res- 
onance frequency of a single isolated meniscus. Super- 
hydrophobic surfaces with particularly large contact an- 
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FIG. 5: Amplitude and phase of an array of micromenisci 
with hexagonal pattern, a — 15/xm and R — 3fim a) at the 
fundamental collective mode (159kHz), and b) at the second 
collective mode (290kHz). Note that the color scale for the 
phase is narrower in a) as compared to b). 

gles or large slip length are expected to show the low- 
est resonance frequencies. Optical diffraction has proven 
an accurate tool to study superhydrophobic surfaces and 
their nanoscopic hydrodynamics. It remains a challenge 
to extend the optical diffraction technique to the accurate 
study of the microscopic shape of the menisci, presum- 
ably by measuring also the angular dependence of the 
diffraction pattern and subsequently solving the inverse 
diffraction problem, following upcoming ideas in theoret- 
ical diffractive optics. 
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